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Abstract 

Quenched reduction is revisited from the modern viewpoint of field- 
orbifolding. Fermions are included and it is shown how the old problem 
of preserving anomalies and field topology after reduction is solved with 



D . 

(— I ^ the help of the overlap construction. 

^ : 1 Introduction 

b 

' The techniques producing numerical values for physical quantities affected by 

the strong interactions have continuously evolved over the last twenty years. 
They have become quite sophisticated and, coupled with major increases in 
affordable computing power, generate steady, incremental progress However, 
if the sole activity in lattice field theory were that of adding digits of higher 
accuracy to some basic quantities one might get the impression that this subject, 
while useful, is rather dull. Luckily, this is not the case: there are lattice 
activities which do take risks and, at times, produce a breakthrough (e.g. lattice 
chiral fermions ||]), or, show that lattice field theory is good for other things 
than QCD (e.g. Higgs triviality bounds [^). There are people who venture 
into supersymmetric territory tackle "insoluble problems" like QCD with 
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a chemical potential Q , use lattice techniques to check scenarios for explaining 
the baryon number asymmetry ||^, and are inventive in many other ways. In 
turn, lattice field theory constructions enter mainstream particle physics [Q 
and even quantum gravity |^] , providing novel "UV-completions" . 

In this talk I shall present a project whose long term objective is to solve 4D 
planar QCD. The project is in its initial stages and, unlike more typical lattice 
field theory projects, is far from being guaranteed to succeed. This project is a 
collaborative effort with J. Kiskis and R. Narayanan At present, the focus 
is narrowed down to the calculation of the Q^-dependent meson-meson correla- 
tion functions in 4D. In 2D, where the planar limit is solvable (as shown by 't 
Hooft jl^), this correlation function is of fundamental importance. Although 't 
Hooft's solution was mainly analytical, to actually get numbers also a numeri- 
cal procedure is required, even in two dimensions. There is little chance for a 
similar analytical breakthrough in 4D. But, we hope that a numerical approach 
- this time containing a stochastic element - will be possible using old ideas of 
large Nc reduction and new lattice fermions. 

't Hooft's planar limit in 4D requires three separate limits, one more than 
in ordinary QCD. One needs to take an ultraviolet cutoff. A, (needed to define 
the theory beyond its perturbative series) to infinity, an infrared cutoff (typically 
the volume of an Euclidean torus), V, to infinity and the number of colors, 
at g'^Nc = A fixed, also to infinity (A is to be traded for a physical scale). The 
natural order of limits is as follows: first V —^ oo, next A — )■ cxd and last Nc —^ oo. 
One knows how to exchange the first two limits without altering the result even 
when the quarks are massless. Less is known about taking Nc ^ oo earlier than 
last, but we shall do that. Otherwise, we are relegated to solving not only QCD 
first, but also QCD-like theories for different numbers of colors, and then the 
relevance of getting the planar limit numerically becomes mostly academic. The 
result is nevertheless important for our project because such investigations 
tell us, quantitatively, how far Nc = oo is from Nc = 3. For us, "solving" planar 
QCD means that a shortcut to the Nc = oo limit was used, and the method 
finds a numerical solution at Nc = oo with considerable more ease than would 
be possible in real QCD. This ancient hope sounds more realistic today than 
it did five years ago also thanks to progress made by string theorists [Q. To 
be sure, in general, the interchange of the limits Nc — > oo with A ^ oo is non- 
trivial iQ. For example, the theory may contain a set of degrees of freedom 
whose number docs not increase with Nc, but whose impact on the ultraviolet 
behavior at any finite Nc is large. Nor is the interchange of the F — > oo and 
Nc oo limits a trivial matter in a theory with massless elementary fields |Tq| . 



2 Large Nc reduction 

We start on a lattice with a "bare" coupling (3 — and at a finite volume 

V. There is very good evidence that the limit iVc — > oo, with b = fixed, 
exists. It exists for the entire range of positive values of b. At small values of 
b (strong coupling) large Nc reduction jl^ tells us that the V oo limit is 
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trivial: there is no F-dependence left at leading order in N^, one can set V — 1 
and lattice momentum space is generated from the eigenvalues of the gauge field 
operators. Complications arise when one tries to extend this to weak coupling 
(large b) but there are some options available Q. This is not dissimilar 

to the situation in models that admit a weakly coupled dual string description 
at large A. 

In their work, Eguchi and Kawai proved that the Schwinger-Dyson equa- 
tions for single-trace operators in the reduced and full models were the same at 
infinite Nc- Later it was realized that this only established reduction for 
small enough b, and that to extend it to large b an additional ingredient, for 
example, "quenching" was needed. 

In modern language, large Nc reduction works by "field orbifolding" |p^ . 
This terminology comes from continuum field theories (actually these theories 
are the low energy sector of some string theories and the term originates from 
a method of constructing special backgrounds for the latter but the de- 

pendence on continuous space-time plays little role in the basic combinatorial 
identity. In our lattice application there is no explicit continuum. Field orbifold- 
ing can produce pairs of theories that admit 't Hooft expansions and, at least at 
the level of bare diagrams, are simply related by possible rescalings of the cou- 
pling constant. The fields in one member of the pair are obtained by subjecting 
the fields in the other member to a set of constraints. These constraints reduce 
the symmetry group from U{Nc) x G to ®iU{Nl) x G'. In the context of large 
Nc reduction one starts with a theory that has F = 1, a single site lattice, and 
one ends up with a theory with a larger lattice (here chosen as a four dimen- 
sional symmetrical torus), V = L^. The original theory has a U{Nc) x U{1)^ 
and the new theory has a symmetry [(8>i=i,L<if/(-/Vc/L'*)] x U{iy x Z{Ly. (The 
N^ — Nc/L'^ and are assumed to be integers.) The two theories have the same 
b and, in the second, the factor (8)4=1,^4 [/(A^^/^'*) is a lattice gauge symmetry 
group. The original theory is the reduced model and the main assumption of 
our project is the hope that the factor isn't really needed to get to the planar 
limit. 

Let me first illustrate the main trick in a simpler model, where the field 
integration measure, excepting the action term, is flat and unrestricted. This 
makes it easy to implement the constraint on the original fields by a "decoration" 
of the original planar diagrams. The "decoration" is a discrete gauge- field theory 
whose gauge link variables need to be summed over, living on the planar diagram 
of the original theory. 

The model is known as the "Weingarten model" |2^, and it defines a gen- 
erating functional for random surfaces on a hypercubic lattice - a rudimentary 
string theory. The partition function of its reduced version ]2l| ] is given by 

^^w(&)=n n 

Although the integral docs not converge, the formula is meaningful as a compact 
definition of an infinite series in b. Similarly, the large N limit of F{h) = 



■V trA^A^A'' A'< -Y^ trA^A^ 



(1) 
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log(^j^^(&)) is also meaningful as a scries in b. It is given there by the sum 
of all planar |<&|^-type Feynman diagrams with a factor of b attached to each 
vertex. 

The reduced model has a symmetry SU{N) x U{lY, where the SU{N) acts 
by conjugation simultaneously on all of the matrices and the [/(I) factors 
independently affect the phases of each A^ matrix. There also is a symmetry 
group made out of elements consisting of permutations of the /i indices 
and conjugations. This "crystallographic group" plays no role in our particular 
application here, and hence will be ignored (it also was ignored a few para- 
graphs earlier). The constrained, or. equivalently, space-time-extended model, 
is defined on a hypercubic toroidal lattice containing L"^ sites. Its partition 
functions is: 

z^nwib) = n n n 

^=1x^=0 l<i,j<N 
E^^. trA^{=o)A^{^+^,)Al{x+t^)Alix)-J2^ Ex ^ (2) 

The symmetry of the extended model is ['E)xU{N')] x U{iy x Z{L)'^. The first 
factor is a local lattice gauge symmetry, the second is a global symmetry and 
the third a space-time symmetry. Again we ignore the crystallographic group. 
Actually, the phases associated with the t/(l)'s go only from zero to ^ - larger 
phases can be reduced by lattice gauge transformations. Again the integral does 
not converge, but the formula is meaningful as a compact definition of an infinite 
series in b'. Similarly, the large N' limit of F'{b') = j^d(]^,y^ log(-^CDw(^0) is 
meaningful as a series in b'. It is again given by the sum of all planar |$|*- 
type Feynman diagrams with a factor of b' attached to each vertex, but, with 
propagators that are less trivial than before. 

The extended model can be viewed as a constrained version of the reduced 
model if we pick N = N'L"^, b = b', write the 1 < i < N indices as double 
indices {i',x) with l<i'<7V',0<x^<L — 1 and restrict the [Afj.](ii^x),(j',v) 
to be zero when y ^ x + jj,. When y = x + fi there are no additional restrictions. 
If all j4^'s in the reduced model of eq. (1) obey this constraint, the portion 
of the trace referring to the x, y components of the indices can be written out 
explicitly and one obtains the extended model of equation (2). An important 
observation is that the constraint can be viewed as enforcing invariance under 
a subgroup of the symmetry group of the reduced model: 

Af, = %{Afj)=e-^^'^''g-^A^g^, u=l,..,d (3) 

[9A(e.x).,(3',v) = ^i',j'^^,v(^^^" (4) 

The set of operations % generates a Z{LY subgroup of the symmetry group 
of the reduced model. The elements of the group are labeled by [s] where 
[s] = (si,S2,...,Sd) and = 0,l,...,i- 1. 

TW(A^) = e^^*- (<?W)^ ^M^'^l, 9^'^ = n at-- (5) 

i/=i 



/ 



dA^Jx)dAf*{x) 
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Note that the action of T^"! on differs from its action on A^, for ^ v by 
an overall phase factor. This Z{L)'^ intersects both original symmetry group 
factors, SU{N) and U{1)'^. Invariance under the action of the four generators 

is equivalent to invariance under the entire Z(L)'^ group. 

Since we integrate freely on the entries of the matrices A^, we can redefine 
the constrained model by keeping full matrices in the quadratic term (where 
the constrained and unconstrained entries of the j4^-matrices do not mix) and 
only matrices satisfying the constraint in the quartic term. The constraint can 
be put in effect by first expanding in the exponential of the quartic factor, and 
then replacing each occurrence of a constrained ^4^, independently, by T[*1(A^) 
with an unconstrained A/j,, and averaging each of the many T^^^ factors over the 
group Z{LY. Only the parts of that are invariant, making up the matrix 
A^j, constr, survive this group averaging, while the other entries are set to zero. 
We now leave this group averaging to the end of the calculation, and, instead, 
do the Gaussian integrals first, producing Feynman diagrams. The propagators 
are now given by 

< T[*1(A,),,TW(4).. >= e^(*''-«'')5!r*'5£-^' V (6) 

T[*1 and T^*! appear only in this propagator. Moreover, as we see, there is only 
dependence on [i — s] . Hence, the independent averages over \t\ and [s] can be 
replaced by a single average over the combination — s] . It is only this extra 
group averaging procedure that distinguishes the two models. 

The sum over the i, j, ... indices produces a factor of for each index loop 
in the double line notation for propagators in the unconstrained model. In the 
extended model, because of the constraints, these index sums produce a trace 
of the product of the various g'l*'! factors encountered round the index loop. The 
gl^'J matrices also form a group, Z[LY, but this group is distinct from the Z[LY 
made out of the T^^l operations. Thus, for every index loop we need to calculate 
the trace of a that is the parallel transporter round the loop. 

From their definition we see that the ^fl**! obey: 

trg^'^ = Ar5[,],[o] (7) 

where the delta function forces the group element to be unity. This has the 
important consequence that the phase factors that appear in the definition of 
the T[*1 operations cancel out any time the contribution from the index loops 
is non-vanishing. (Note that the phase factors are associated with propagators, 
not index lines.) Thus, we can ignore the phase factors and take the propagators 
as given by 

<TW(4),,TW(4),,>-5!Mr' (8) 

Group averaging amounts to averaging over all L'^ values of the [u] labels. 

More generally, this induces us to think about a situation in which we assign 
to each propagator in the planar diagram a group element g € G where G is some 
finite group, g is represented by an A'' x A'' matrix Dij{g) which replaces one 
of the two index lines making up the propagator. Before the "decoration" , this 
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index line was a Kronecker delta. The other index line of the same propagator is 
decorated by Dki{g~^)- The group element associated with each propagator is 
independently averaged over all of G and the decoration introduces a weighting 
factor given by the product over all index loops (faces) of traces of products of 
D{g) matrices round the loop. The summations over G are now easy to perform 
- a standard exercise in lattice strong coupling expansion, this time defined on 
a random two dimensional lattice |22| . 

For any h G G, trD{h) is a class function and hence expandable as: 

trD{h) ^^CrdrXr{h) Cr = ^Xr{h^^)^trD{h) (9) 

The averaging over the g^s associated to the propagators produce a total factor 
given by 

Y^cidi = Y.-l^''di (10) 

r r 

where / is the number of faces and dr the dimension of the r-irreducible rep- 
resentation of G. In the above equation we make use of the planarity of the 
diagram. If e is the number of edges in the diagram and v the number of vertices, 
we have: 

f-e + v = 2 e = 2v f^2 + v (11) 

In order to be able to absorb the decoration factor into a rescaling of b' we 
need Cr — K for all r. From (9), K — ^j^, where d{G) is the number of 
elements in G. d{r) is the dimensionality of the irreducible representation r, 
and d{G) = Ylir ^r- The contribution of the diagram is (^jr^ J2r ^ is an 

integer and we see that the representation of G furnished by D{h) contains K 
copies of the regular representation. This is a consequence of the requirement to 
have an overall b''" factor. It is easy to check that in our case, where G = Z{L)'^ 
(with elements g^'^^, this condition on the representation holds. 

Turning now to our special case, we set d{G) = L"^, N' ~ and b' = b. 

The weight of the above diagram becomes N^b^-j^- From the definition of the 
free energy of the reduced and extended model we now conclude that 

F'ib) = ^ log(Z^Dw(^')) = m (12) 

The equality of the free energies is the main result of reduction. Relations 
between other correlation functions can be obtained by adding more terms to the 
action and taking derivatives with respect to the new, infinitesimal, couplings. 
The above result was obtained a long time ago by direct inspection of the 
series and by mapping the Schwinger-Dyson equations of single-trace operators 
from one model to the other. 

In its extended version, the Weingarten model will turn into a lattice gauge 
theory if we constrain the A^{x) matrices to be simple- unitary and change the 
integration measure to the SU{N) group Haar measure. The relation between 
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the reduced and extended versions holds even after the unitarity constraint has 
been imposed on both models. This was proven a long time ago, first by the 
Schwinger-Dyson equation method p6[ | and, subsequently, by inspection of lat- 
tice strong coupling diagrams [ p3| and also using stochastic quantization p^ . 
A proof by field orbifolding would require a generalization of the analysis above, 
but since the result is known to hold, there is little doubt that such a generaliza- 
tion exists. Field orbifolding admits generalizations involving nonabelian global 
symmetries, and it would be interesting to see whether previously unknown 
results in planar lattice gauge theory could be obtained from this approach. 

3 Prototype model and its extended variants 

The unconstrained, reduced, gauge model consists of 4 matrices G SU{Nc), 4 
unitary matrices S SU{Nf) and two rectangular Grassmann matrices ^, 5", 
where one index is in the fundamental representation of SU (Nc) (color), and the 
other index in the fundamental representation of SU{Nf) (flavor). The fermion 
matrices also have a Dirac index, which we often suppress. The model is defined 
by: 

5 = E ^^-^^^i +f^'T. + Tr [^D{U, V) (*)] 

c;,^[uM ci,^[vM (13) 

The linear operator D{U, V) is a reduced version of a lattice Dirac operator, to 
be introduced later on. 

The constrained associated model lives on an L"^ lattice . The constraints 
mean that the matrices C/^^, V^, 'I', 4' are expressible in terms of smaller ma- 
trices Ufj,{x)^Vfj,{x),ip{x)^ip{x). The dimensions are ric instead of Nc and n/ 
instead of Nf where ricj = Ncj/L'^- {ficj, Ncj, L are all integers.) 

\-^tJ'\ix),{jv) ~ Wfi{^)]ij^y,x+tJ.7 i^P-]{ax),{by) ~ [^fi{^)]abSy,x+tJ. 

[^]i^x),iay) ^ i^'(^)]^-Sxy, = [V^(a;)]-'5,. (14) 

For small lattice gauge coupling the constrained and unconstrained models will 
be equivalent in the limit Ncj oo with ^ held fixed and a simple D{U, V). 
But, for larger coupling the equivalence breaks down and we shall adopt here one 
of the proposals for keeping the equivalence up even at large lattice couplings, 
namely quenching fl^ . 

4 Quenching 

Quenching amounts to freezing the eigenvalues of ?7^, to sets that are uni- 
formly distributed round the circle. Other than that the unitary matrices are 
free. One also takes the rescaled coupling, j^, to infinity so as to freeze out the 
flavor gauge group p6[ | . The C/^ eigenvalues play the role of lattice momenta in 
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boson propagator lines, while the eigenvalues enter only into the fermionic 
lines. The operator D{U, V) breaks up into Nf blocks of size iNc x ANc each. 
This is the crucial point of our approach: The numerical price of handling 
fermions has a chance to simplify tremendously if indeed Nc does not need to 
be anywhere nearly as big as the parameter ricL'^ typical to numerical QCD 
simulations and, nevertheless, the planar limit stays well approximated. 



5 Overlap fermions 

We now define D for one of the flavor blocks. So long we are focusing only 
on fermions we can absorb the flavor phase into the color matrices C/^i (whose 
eigenvalues are quenched as before). Define the unitary operators X^: 



X, = '-^U, + i±^[/t (15) 



The Wilson-Dirac matrix Dw is a given by 

Dw = m + Y,i^-X^) (16) 

We choose the parameter m to be about —1 (it must be between —2 and 0). It 
is easy to see that Hw = 75-Dvy is hermitian. 

The overlap construction says that the chiral determinant line bundle is 
represented on the lattice by the inner product of two states, (+|— ). Here (+| 
is the eigen-bra of a system of noninteracting second quantized fermions in a 
single body potential given by Hw with the mass sign switched, while |— ) is 
the eigen-ket of a system of noninteracting second quantized fermions in a single 
body potential given by Hw with the mass chosen as above. One can change 
the lattice definition without affecting the continuum limit so that the mass in 
Hw is taken to be +00 for the (+|-state. In that case the Hw for the (+|-state 
can be replaced by 75 |2| . 

The overlap (+|— ) can be written in terms of single particle wave- functions 
as the determinant of a matrix A defined below j^, Choose the chiral 

representation for the 7^ matrices, so that 

7.^('°,) (IT) 



^ -1 ^ 

Let $ be a simple unitary matrix of eigenvectors of Hw, ordered so that 

sign(i?M.)$ = $ J 0^ ) , '^-{c d) (1^) 

Then 

(+|-)=detA (19) 
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For a vector-like theory we need | det(A)p, which, unhke det A, is a function, 
not a section of a bundle. Simple linear algebra tells us that 



However, 



Hence, 



det 



I det (A) |2 = det (A) det 
1 + 75sign(i7n/ 



A -B 
-C D 



det 



1 + 75sign(77, 



w) 



(20) 
(21) 

(22) 



We discover that chiral symmetry can be preserved on the lattice if we use the 
overlap Dirac operator M: 



Dn = 



1 + 75signi/vv 



(23) 



The operator D is constructed out of Do, with the addition of a small mass 
term for the quarks. 



6 Topology and Anomalies 

Anomalies are expressed by particular Berry-phase factors innate to the |— ) 



state of the overlap |28 . They obstruct smooth gauge invariant sections in 
cases where anomalies do not cancel. 

Topology is expressed by H-^j^ having unequal numbers of positive and nega- 
tive eigenstates (above, where A is taken as a square matrix, we tacitly assumed 
these two numbers to be equal but this need not be so, in which case what we 
called "det A" is identically zero) . One can easily exhibit explicit sets of 
matrices for which there is non-zero topology. 

Problems arise if Hw has an eigenvalue that is exactly zero. But this will 
not happen as the continuum limit is approached because one can prove the 
following inequality poj : 

\\\U,, U^W < e =, MS^) > 1 _ |i + „,| _ Sie) (24) 

where (5(e) goes to zero as e goes to zero. 

At infinite Nc, for large enough lattice couplings and a standard lattice 
action, the commutators indeed become bounded j^, as required. 

7 Some history 

The first treatment of fcrmions in reduced models was carried out by Gross 
and Kitazawa [pO| who invented the momentum "force-feeding" prescription 



9 



that was derived above as a remnant of the frozen-out flavor gauge sector . 
These authors also invented a version of reduction that used hermitian matrices, 
Afj^, rather than the unitary link matrices one has on a lattice |3^. The Gross- 
Kitazawa model can be obtained from the lattice model by approximating w 
1 + lA^. However, as also noted by Gross and Kitazawa, this approximation 
does away with anomalies and topology. For many years we did not known how 
to cure this problem on the lattice, where the link matrices are unitary. The 
problem was the same as the notorious problem of lattice chirality. As we have 
seen, recent progress on the chiral fermion problem |^ makes it now possible 
to preserve all the necessary fermionic properties also in reduced lattice models 
removing the main obstacle of principle to making nontrivial use of reduction 
in numerical simulations. 
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